Dynamical quantum phase transition of a two-component Bose-Einstein condensate in 

an optical lattice 
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We study dynamics of a two-component Bose-Einstein condensate where the two components are 
coupled via an optical lattice. In particular, we focus on the dynamics as one drives the system 
through a critical point of a first order phase transition characterized by a jump in the internal 
populations. Solving the time-dependent Gross-Pitaevskii equation, we analyze; breakdown of adi- 
abaticity, impact of non-linear atom-atom scattering, and the role of a harmonic trapping potential. 
Our findings demonstrate that the phase transition is resilient to both contact interaction between 
atoms and external trapping confinement. 

PACS numbers: 03.75.Mn,64.60.Ht,64.70.Tg,67.85.Hj 



I. INTRODUCTION 



Since the pioneering experiments on Bose-Einsten con- 
densates (BEC) [1], the field of ultracold atomic gases has 
seen a tremendous development Nowadays, prepa- 
ration and manipulation of BECs is a standard proce- 
dure providing a versatile testbed for the study of var- 
ious quantum effects. Spinor condensates, where inter- 
nal atomic Zeeman levels play an important role for the 
dynamics, have been experimentally studied by numer- 
ous groups More recent experiments on spinor con- 
densates include; coherent transport spin-mixing 5|, 
inherent spin tunneling Q, and symmetry breaking 7|. 
Placing the condensate in an optical lattice formed by 
counter propagating laser beams greatly affects its char- 
acteristics and exciting phenomena, e .g. Bloch oscilla- 
tions P , gap and colliding solitons d, ^ , self-trapping 



ll| , superfluid instability [I4I , vortices [13| , Anderson lo- 



calization [3], and for more strongly correlated systems 
the superfluid-Mott insulator phase transition (PT) [3] , 
arise. 

Most theoretical works studying PTs consider systems 
at strict thermal equilibrium. However, a more appropri- 
ate picture capturing the underlying physics of real ex- 
periments is often encountered by a dynamical approach. 
Indeed, the great experimental progress seen in the recent 
past on correlated ultracold atomic systems [TgI, calls 
for a deeper theoretical understanding of non-equilibrium 
PTs. Much of the theoretical works are devoted to dy- 
namics across a critical point in various spin models 
[S m, d, S Also the dynamics of phase tran- 

sitions in lattice many-body systems such as the Bose- 
Hubbard model have been considered [13, H^. To de- 
rive an estimate of the adiabatic breakdown, the Landau- 
Zener formula, giving a measure of transition probabili- 
ties across an avoided level cros sing (2 4 1 , has been utilized 
in most of the above works flst 1^ 20l . [2^ . Moreover, 
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in Refs. [17|, [IJ, OJ, [22, l23| the Kibble-Zurek mechanism 
of quench induced transitions was employed. 

In this paper we consider a spinor condensate in an 
optical lattice. Apart from providing an effective peri- 
odic potential, in this model the lattice also induces a 
coupling between the internal atomic states. The cor- 
responding many-body model, without atom-atom inter- 
actions, was first introduced in Ref. [2^ . and the occur- 
rence of a first order PT between different internal states 
was demonstrated. The system was further investigated 
in by taking interaction between the atoms into ac- 
count. Both of these works consider thermodynamical 
equilibrium. Since the single particle ground state of the 
present model is not gapped, even small time-dependent 
perturbations of the Hamiltonian might well cause non- 
adiabatic excitations. Driven through the critical point, 
such excitations might wash out the signatures of the 
PT. One goal of the present paper is to study the im- 
portance of these non-adiabatic contributions. Moreover, 
the model automatically takes into account the effects de- 
riving from atom-atom interactions, and the state chang- 
ing collisions they bring about. This aspect was left as 
an open question in ^25]. In addition, the sensitivity to 
external harmonic trapping is also considered. 

The outline of the paper is as follows. In the next 
section, we present the one dimensional single particle 
Hamiltonian. After first discussing some general prop- 
erties in Subsec. Ill Ai we numerically diagonalize the 
Hamiltonian and give the spectrum in Subsec. IIIBI Us- 
ing the spectrum, we demonstrate the equilibrium PT 
in Subsec. Ill CI The following Sec. IIIII is devoted to our 
main results; the dynamics when the system is driven 
through the critical point. The numerical method for 
solving the Gross-Pitaevskii equations is outlined, and 
first the thermal equilibrium PT is studied by means of 
the partial-state fidelity susceptibility. In Subsec. IIII AI 
we calculate the order parameter for different dimensions, 
whereby the following Subsection considers the situation 
of external harmonic trapping. Breakdown of adiabatic- 
ity is analyzed in Subsec. IIII CI Finally we conclude in 
SecHVl 
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II. ONE DIMENSIONAL IDEAL GAS MODEL 
SYSTEM 

The appearance and nature of the PT are easily ex- 
tracted from the structure of the spectrum of the one 
dimensional single particle Hamiltonian. The effects of 
particle interactions, higher dimensional systems, and ex- 
plicit time-dependence of system parameters will be con- 
sidered in Sec. IIIII 



A. Ideal gas Hamiltonian 

The properties of an ideal gas are captured in the cor- 
responding single particle system. We consider a three- 
level A atom dipole coupled to two light fields as illus- 
trated in Fig. [T] The 1^3 transition is driven via a pe- 
riodic standing wave field, while an "external" spatially 
independent laser gives rise to the 2 <-> 3 transition. This 
assumes the perpendicular laser to have a mode-waist 
larger than the extent of the condensate. The effective 
atom-field couplings are A and Q, the wave number fc, 
and the atom-field detunings Si and 62 respectively. In 
the case of large detunings, the excited atomic state |3) 
may be adiabatically eliminated yielding an effective 2x2 
Hamiltonian [23 



Hi 



2m 



hA 



-az — hUi cos(2fca;)CTii -|- hU cos{kx)ax, 

(1) 

where A = — ^2! — ^'^/S2 — A^/25i is an effective 
detuning taking into account for the constant Stark shifts 
of states 1 and 2, Ui = X^/2Si, U = Xn{l/26i + I/2S2), 
and the tr-operators are the Pauli matrices; (t^ — \2){2\ — 
= |l)(2|-t-|2}(l|, and<Tii = |1)(1|. Respectively, 
p, X, and m are atomic momentum, position and mass. 
Note that the Ui term originates from the Stark shift of 
the |1} state. As it affects only the |1) state and not the 
2) state, it shifts the detuning resonance Ac = ^ Ac 7^ 
0. Apart from altering the detuning resonance, it also 
renders a light induced potential for the state The U 
term, on the other hand, provides an effective potential 
for both internal states. In the limiting cases, when one 
of the parameters A, Ui and U is dominating the other 
two, it is legitimate to assign a particular potential to 
the internal atomic states. In general, however, the full 
coupled system has to be considered as an entity [1^. 
For later convenience, scaled quantities are introduced 

,-2.2 

by the characteristic energy and length scales Er — 
and k~^, such that 



U,,'^. a = f. (2) 

rjr i^r 



In the matrix-form, the scaled Hamiltonian reads 

A 
2" 



Hi 



Ucos{x) 



U cos{x) 



A 



Ui cos(2a;) 



(a) Atomic condensate 
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Optical lattice 



(b) Atom-laser configuration 
13) 




External laser 



FIG. 1: Schematic configuration of system setup (a) and 
atom-laser configuration (b). An optical lattice drives the 
|1) ^ |3) atomic transition and the "external" laser the 
2) ^ |3) transition. 



where we have |1) = 



and |2) 



From the 



fact that cos(£)|(7} — {\q -|- 1) -|- |g — l))/2, where p\q) — 
q\q), one finds that a swapping of the internal states 1 
and 2 is accompanied by a shift of momentum by one 
unit in either positive or negative direction. Moreover, a 
given momentum eigenstate \q) is only coupled to other 
momentum states \q + rf), where rj is any integer. In 
particular, the states (* = 1:2) can be divided 

into two sets 



\^ri{q)) = 

\^vi<l)) = 



rj even 
77 odd 

T] even 
rj odd. 



(4) 



which are not coupled by the Hamiltonian ([3]); 
{ifri'{q')\Hi\4>ri{q)) = 0. Thus, the Hamiltonian may 



be written as i?i — H^ 



Ha,, where the two sub- 



Hamiltonians operate on states spanned by {\^ri[q))} or 
{\4>ri{q))} respectively. Here, the quasi momentum q re- 
sides in the first Brillouin zone. 

As a periodic operator. Hi commutes with T = e^*'^, 
where ^ = 27r is the scaled period. It is readily shown [2^ 
that the operator 



1 = ^, 



(5) 



defines another symmetry, describing a half-period boost 
of the position and an atomic inversion. An outcome of 
this 1/2 symmetry is that the Brillouin zone extends be- 
yond (—tt/1. tt/1] as imposed by the Z-periodicity of the 
system [25l . |28|. Thereby, the first Brillouin zone is in- 
stead defined for quasi momentum q e ( — 1, 1]. 



B. Ideal gas energy spectrum 

As was already pointed out, the character of the PT is 
directly seen from the energy spectrum. Due to the block 
(3) diagonal form of the Hamiltonian, Hi = H^ (g) H^, the 
dispersions of Hu> and H^ can be calculated separately. 
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In the figures, we will differentiate between dispersions of 
and i?^ by using dashed or solid curves, respectively. 

The internal two-level structure of the system brings 
about anomalous dispersions that possess several minima 
[2^ [2^ . The lowest energy of each band Ei, (q) , where 
1/ = 1, 2, 3, ... is the band index and q e (—1,1] the 
quasi momentum, is twofold degenerate. However, the 
degeneracy is due to the fact that the individual energy 
bands El'^\q) — E\f\q') for some q and q' . The only 
exception is at exact resonance A = Ac, where the bands 
e\^'^ (q) and eI'^'' (q) become identical. 

The Hamiltonian ([3]) is readily diagonalized by trun- 
cating the number of basis states If Ui ^ U , atoms 
in the |1) state will feel a strong periodic potential, while 
2) atoms experience only a weak potential. Thus, the 
spectrum is in this case typically made up of a fairly 
flat band mixed together with "parabolic" energy bands. 
When U dominates Ui , the situation is different and mix- 
ing of flat and "parabolic" bands does not exist. Exam- 
ples of the first three energy dispersions are displayed in 
Fig.d In (a)-(d) U > Ui, while in (e)-(h) Ui > U. In 
the last three plots, almost flat dispersions are found due 
to the large Ui, and especially in (f) there are "parabolic" 
bands with lower energy than the flat one. 

C. First order phase transition 

At zero temperature, the condensate atoms will oc- 
cupy the ground state. Even though the ground state 
is degenerate, mainly one of the degenerate states will 
be populated. In an experiment, the atoms will typi- 
cally be initialized in either jl) or [2) and then cooled 
down to its ground state. Hence, the obtained initial 
ground state will predominantly populate either of the 
two energy bands E['^\q) or E['^\q). An equal popu- 
lation balance between the degenerate ground states is 
very unlikely, and indeed, in our numerical simulations 
of interacting systems utilizing the imaginary time prop- 
agation method we never encounter such balanced situa- 
tions. 

Comparing Fig. [5] (a) and (d), it is seen that the dis- 
persions E'lf\q) and Ei'^\q) have been roughly swapped 
when A ^ —A. Thus, the minimum of E'"^\q) or 
E^\q) are changed between the center and edges of 
the Brillouin zone as A is tuned through the resonance 
A = Ac. In the example of Fig. [5] (a)-(d), this resonance 
condition is Ac « since Ui <^ 1. When Ui is large, 
however, the Ac is shifted away from zero, but nonethe- 
less, between A < Ac and A > Ac the minima is either 
at the Brillouin center or at its edges. As the quasi mo- 
mentum corresponding to the energy minima of either 
E'f\q) or E^\q) changes abruptly when A is tuned 
across the resonance, one encounters a first order PT. 

As discussed in the previous subsection, the internal 
(|1) and |2)) and external (momentum) degrees of free- 
dom cannot be changed independently via the Hamil- 



tonian ([3]). Absorption or emittence of a lattice photon 
always induces a momentum "kick" . Consequently, shift- 
ing the momentum of the ground state between the cen- 
ter and the edges of the Brillouin zone is only possible 
by swapping the internal states |1) and |2). Thus, the 
collective inversion for atoms belonging to either of the 
E^;^\q) or E[^\q) bands 

ji^^-j;^T.('^i'^)^ 7 = '^,</', (6) 

where al is the ith atom's z Pauli matrix, should be dis- 
continuous at Ac. Here, Na is the number of condensate 
atoms in band E['^\q), and we note that since we con- 
sider only the ground state, (o'i*'') = (iri"''') for any i and 

j. The full inversion jz = ji^^ + jl'^'', however, depends 
on the number of atoms residing in each band. 

Figure [3] displays two examples of jz when all atoms 
belong to the £'|'^^(g) band. In (a), the inversion is shown 
for the parameters corresponding to Figs. [2] (a)-(d), and 
the gap of the inversion is found to be well resolved. In 
(b), however, corresponding to Figs. [H (e)-(h), reveals a 
small gap in the inversion. In this case, the PT occurs for 
a detuning Ac ~ 3.8 which is somewhere between Fig. [2] 
(g) and (h) where the dispersions are fairly flat. We note 
that apart from shifting the critical resonance Ac, large 
JJi's also results in asymmetric jl'^'''s with respect to Ac. 



III. DYNAMICAL PHASE TRANSITION OF A 
WEAKLY INTERACTING BOSE-EINSTEIN 
CONDENSATE 

Section |TT] demonstrated that a first order PT is en- 
countered by driving the system of non-interacting atoms 
through a critical detuning Ac. Adding interaction be- 
tween the atoms will introduce coupling between the two 
sets of basis states (|3]). This effect is, however, relatively 
small f26l]. On the other hand, very strong atom-atom 
scattering amplitudes can have other important conse- 
quences, for example it may drastically change the sin- 
gle particle dispersion curves [s^]. Nonetheless, it was 
argued in Refs. [25l . [26j . that the PT should be stable 
even against strong interaction between the particles. In 
this section we will show, on a mean-field basis, that 
this is indeed true, and that the PT is also insensitive to 
an external harmonic trapping potential. Moreover, the 
dynamics of the condensate driven through the critical 
point will be thoroughly addressed. 



A. Sweep through the critical point 

The spectra E^^\q = 0, A) and E^^\q = ±1, A), and 
the spectra E^^'^ {q = 0, A) and E^f' {q = ±1, A) become 
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FIG. 2: Examples of the three lowest energy bands of Hamiltonian (|3]). In (a)-(d) we have U > Ui {U = 0.1 and Ui = 0.05), 
while in (e)-(h) U < Ui {U = 0.2 and Ui =8). The multiple number of local minima of the lowest band is clear. Dashed lines 
are the energy dispersions El^\q) of the sub-Hamiltonian H^p, and solid lines are E[f\q) corresponding to H^. 



(a) 1 







U=l 
Uj=0.5 







-2-1 1 2 




3.5 4 4.5 

A 

FIG. 3: The inversion (|6} as function of detuning A. In the 
upper plot, t/i <C 1 and Ac ~ 0, while in the lower plot f/i = 8 
giving Ac « 3.8. 



degenerate and cross at A = Such crossing of ener- 
gies is characteristic for a PT [3l|. The crossing of ener- 
gies in the present model is envisaged in Fig. 3] displaying 
the energies E^\q,/S.) as function of A for q = —n/N, 
where A'' = 10 and n = 0, 1, ... N. Letting N ^ oo, the 
spacings between the lines in Fig. [3] vanish. The two fig- 
ures use the same parameters as in Fig. [31 and we note 
that when Ui is large, Ac is shifted away from zero as 
was already seen in Fig. [3] Further, in (a) corresponding 
to Fig. [3] (a) where the discontinuity of the order param- 
eter is more pronounced, the energies are more sparsely 
spaced compared to (b). Indeed, in (b) it is not even clear 
that the energies e['*'^ {q = 0, A) and e[''''^ (g = ±1, A) ac- 
tually cross. However, a closer look around the critical 
point reveals the crossing. 

The purpose of this section is to study the effects on 
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(b) -1.8 




FIG. 4: Examples of the lowest band energies El '{qn, A) for 
q„ — — n/10, n = 0, 1, ... 10. The parameters are the same 
as in Fig. E] [/ = 1 and f/i = 0.5 in (a), and U = 0.2 and 
I7i = 8 in (b). The dotted line and the dot-dashed line are 
the energies E['^\—1, A) and e['^\o, A) respectively, and the 
crossing between these two, e['^\o,A) = — 1, A), gives 

the critical detuning Ac. 

the PT arising from interaction between the particles and 
from a confining harmonic trap. Furthermore, the analy- 
sis will not be restricted to the one dimensional case, but 
we also considers the two dimensional counterparts. One 
interesting aspect is the dynamics in itself as the system is 
driven through the critical point. Another central aspect 
that will be investigated is breakdown of adiabaticity. 
Due to the continuous spectrum and the crossing of the 
energy levels, non-adiabatic excitations of the conden- 
sate are expected and especially around A — A^, where 
the density of states drastically increases. 

In the higher dimensional cases we assume, for simplic- 
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ity, that all lattice fields couple with the same strength 
to the atom and that they share the same periodicity. 
Such a scheme requires multiple excited Zeeman levels 
as well as mutually orthogonal polarizations of the lat- 
tice beams, but after the adiabatic elimination on all of 
them, one is left with an effective two-level Hamiltonian 



where 



A 

2 



U^COs{Xa) 



U'^^cos{xa) — — — C/i^cos^(a 



(7) 



Here a indicates a component of the vector along the 
direction a. 

Moreover, we will take into account for collisions be- 
tween the particles in terms of a non-linear term in the 
Hamiltonian, as well as consider a harmonic trapping po- 
tential. We thus consider the Gross-Pitaevskii equation 

«^«'(x,<) = i/GP*(x,i) 



;2 -f 5't(x).g.5-(x) 



*(x,t), 
(8) 

where t is the scaled time, uj the scaled parameter de- 
termining the trap frequency, and g the scattering am- 
plitudes. In this coupled two-level problem, the wave 
functions are spinors 



^i(x) 



and 



.911 512 

9l2 922 



(9) 



(10) 



Normalizing the wave function as / (ix|^(x)p = 1, the 
scattering amplitudes become gij — A'KhNaij/mVEr^ 
where N is the number of condensate atoms, V the ef- 
fective volume, and the state dependent s-wave scat- 
tering amplitude. 

We will solve Eq. ([8]) numerically, starting in the 
ground state for large detuning and then change A in 
time such that it passes the critical point. In particular 
we chose a linear sweep of the detuning. 



A(0 = Aot, 



(11) 



where Aq sets the sweep velocity. We integrate (|5]) from 
large negative to large positive times t G \—T, T] such 
that the critical point is surely traversed. 

The GP equation ( [8|) is solved using the split operator 
wave packet method j32l |. The initial state is taken as 
the ground state for the given system parameters. In 
particular, the ground state is obtained via the imaginary 
time evolution 

vI/G5(x) = hm , "^"^""'^^ , (12) 
dx|*o(x,r)|2 



^'o(x,t) = e 



*o(x,0) 



(13) 



and (x,0) is the initial wave function ansatz which 
will be taken as a normalized real Gaussian. The sub- 
script indicates the initial Hamiltonian parameters U, 
Ui, g, Aq, and T, as well as the center of the Gaussian 
xq = and its width Ax (which is taken to be much 
larger than the period Z = 27r of the Hamiltonian). This 
numerical procedure gives a good approximation for the 
ground state. Another important point is that since the 
ground state (for a non-interacting system) is degener- 
ate, the imaginary time evolution renders different states 
depending on the ansatz (|13p . In particular, since for the 
state of our assumption all population resides in the |1) 
state and the momentum is centered around the origin, 
the ground state obtained from (jl2p will mainly popu- 
late states spanned by {\Lp,j{q))}. This holds even for 
relatively large g. Again, an imbalanced initial state is 
believed to be in good agreement with true experimental 
realizations; roughly speaking, the imaginary time evo- 
lution represents the cooling of the condensate. 

Before studying the actual dynamic problem, we con- 
sider the FT characterized by the ground states numer- 
ically obtained via imaginary time evolution (|12p . By 
doing so, we get an idea of how sensitive the FT is to 
the non-linear atom-atom interaction term. In particu- 
lar, we evaluate the partial-state fidelity susceptibility, 
which has been demonstrated to give a clear indication 
of various FTs [s^. The advantage of utilizing a fidelity 
measure, originating from the quantum information con- 
cepts, to analyzing FTs is that it is state independent 
and no prior knowledge about order parameters is re- 
quired. The fidelity susceptibility has been shown to con- 
tain sufficient information to reveal the universality class 
of the FTs [sl, even for topological FTs [sl, Kosterhtz- 
Thouless transitions and FTs characterized by an 
underlying continuous level crossing (36j . The idea is sim- 
ple, letting the ground state be parameterized by some 
quantity h (in our case h = A), p{h) — |^'(/i)}(^(/i)| we 
introduce the reduced density operator for the internal 
states paW = TiB[p{h)], where the trace is over ex- 
ternal degrees of freedom. For the present model, it is 
straightforward to derive 



N2i{h) N22{h) 



(14) 



where iVy 



„iiv.iv. Nij = J dxip*{x,h)^pj{x,h) and ipi{x,h) is the 
constituent ground state wave function for the parameter 
h. The partial-state fidelity is defined as 



FA{h,dh) ^ Tr 



\/ pA{h)pA{h + 6h)y/pAih) 



(15) 



In terms of FA{h,5h), the partial-state fidelity suscepti- 
bility reads 

-2\-siFA{h,5h) 



X{h) 



lim 



5h^ 



(16) 
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For the situation at hand, there is a critical hc{= Ac) for 
which the PT occurs. The fidehty measures the overlap 
between the groundstates for different parameters h and 
h + Sh, and whenever he lies between these two values 
the fidelity will drastically drop and the susceptibility in- 
crease. This is displayed in Fig. \5\ The upper plot (a) 
represent the case with U > Ui, while (b) shows examples 
of the reverse situation. For solid lines, the atom-atom 
interaction is zero and for dashed lines it is non-zero. 
The sharp peaks of the susceptibility at he coincide well 
with the critical parameters found in Fig. [3] using the full 
numerical diagonalization of the single particle Hamilto- 
nian. The slight discrepancy is believed to derive from 
the fact that the imaginary time propagation is stopped 
after some finite time t. Interestingly, the shape and lo- 
cation of the peaks of x(^) seems to be insensitive for 
atom-atom interactions, proving that the PT is not re- 
stricted to an ideal gas. 
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FIG. 6: Examples of the collective inversion jz as function of 
time; one dimension (solid lines) and two dimensions (dashed 
lines). In (a) [/ = 1 and Ui = 0.5, while in (b) U = 0.2 and 
Ui = 8. In both plots Ao = 10, lu — 0, gu = (?22 = gi2 = 0. 
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FIG. 5: The partial-state fidelity susceptibility ([T6|l for the 
one dimensional model. In (a) [/ = 1 and Ui = 0.5, while 
in (b) U — 0.2 and Ui — 8. In both plots cu — 0, and gu = 
g22 — gi2 = (solid lines), gu = 322 = IO512 = 10 (dashed 
lines). 



Let US now turn to dynamics. Once the ground state 
is initialized, the wave packet is let to evolve using the 
real time propagator. Since the Hamiltonian is explicitly 
time-dependent, the evolution is non-trivial. Figure [5] 
shows examples of the collective inversion jz = +jz^^ , 
corresponding to the parameters of Fig. [31 Solid lines 
give the one dimensional case and dashed lines the two 
dimensional case. The figure makes clear that the PT 
is resilient to the dimensionality of the system. In fact, 
the change in the inversion is even more pronounced in 
the two dimensional situation. The shift of the critical 
point between the two dimensions seen in (b) can be un- 
derstood from the Hamiltonian ([7]). Since in two dimen- 
sions the sum on the diagonal contains two terms and 
both terms shift the resonance condition, the shifting is 
larger for higher dimensions. 



B. Effect of a trapping potential 

The PT demonstrated in the previous section is an out- 
come of the periodicity of our system and in particular 
the anomalous energy dispersions. The band structure of 
the spectrum will be lost if the condensate is placed in an 
additional external trapping potential. For weak confine- 
ment, however, the spectrum still shows similarities to a 
band-gap spectrum. More precisely, if the characteris- 
tic length scale of the harmonic trap is much larger than 
the lattice periodicity I — 27r, signatures of the PT are 
supposedly still present. On the other hand, very tight 
trapping may wash out any signatures of the PT. 




FIG. 7: The collective inversion jz as function of time for 
the one dimensional case. Solid line {uj = 0.01), dashed line 
{lu = 0.02) and dotted line {uj = 0.05). The rest of the di- 
mensionless parameters are, U = 1, Ui = 0.5, Ao = 5, and 
311 = S22 = gi2 = 0. 

In Fig. [7] we display examples of the collective inversion 
for three different trapping strengths. For these parame- 
ters the three lines are almost identical, and still for the 
largest u the ground state wave function only extends 
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over roughly five lattice periods. By increasing w fur- 
ther one cannot talk about a periodic system, and it is 
found that the momentum wave packet does not undergo 
a sudden momentum shift when driven through the crit- 
ical point. The atom-scattering g = in these examples. 



C. Breakdown of adiabaticity 

As already pointed out, the ground state is not gapped 
in the current system, which should make it sensitive to 
non-adiabatic effects. Adiabaticity is still most likely to 
break down at the critical point where the density of 
states becomes the highest and the instantaneous en- 
ergies change rapidly. In typical models, adiabaticity 
is expected to be violated in the thermodynamic limit 
TV — > oo, since then the curve crossing at the critical 
point becomes infinitely sharp [3l| . Here, we encounter a 
different situations because already at the single particle 
level, the energies at the crossing are degenerate making 
it a non-avoided crossing. Thereby, the constrains for 
adiabaticity can be highly tested already for small atom 
numbers, interpreted as small scattering amplitudes g. 

In this section we analyze the excitations caused by 
non-adiabatic effects. The time evolved wave packet 
reads 



(17) 



where T is the time-ordering operator and ^'gs'(x) is 
given in Eq. (|12p . At any instant of time we can apply 
the imaginary time evolution to find the corresponding 
ground state ^'^^(x) by simply substitute ^(x, t) into 
Eq. (fT^ . We introduce the non-adiabatic relative energy 
increase as 



AE{t) 



E^it)-E^,,> it) 



(t)\ 



(18) 



Aq. The breakdown of adiabaticity is dominantly oc- 
curring around the critical point as expected, and more- 
over, fast sweeps gives a larger breakdown. The case 
of a non-zero scattering is depicted as a dot-dashed line 
in the figure. For this choice of scattering amplitudes 
(511 = 522 = 10(7i2 = 20), the adiabaticity is approxi- 
mately unchanged. However, our numerical results in- 
dicates that around 511 = 1722 = Wgi2 = 30 a drastic 
drop in the adiabaticity appears. Such sudden decrease of 
adiabaticity has also been demonstrated in related mod- 
els M. 
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FIG. 8: The relative energy increase (|18p for different sweep 
velocities (dotted line: Ao = 1, solid line and dot-dashed line: 
Ao — 10, and dashed line: Ao = 50) and scattering ampli- 
tudes (dotted/solid/dashed lines: g = and dot-dashed line: 
ffii ~ <?22 = 10(712 ~ 20). Increasing Aq decreases the adi- 
abaticity as predicted. The solid and dash-dotted lines cor- 
respond to the same sweep velocity but different scattering 
parameters g, and it is noted that the non-linearity does not 
destroy adiabaticity in this parameter regime. However, for 
larger non-linearity we have found great decrease of adiabatic- 
ity. Note that we have scaled time with the sweep velocity 
on the a;-axis. The other dimensionless parameters are as in 
Fig. [3] (a). 



where E^j,{t) and E, (f) (t) are the system energies given 

by the functional E[i!] and respectively. The 

condition dAE{t) /dt ^ 1 is an indication that the evo- 
lution is adiabatic. 

We note that by neglecting the spatial variation of the 
optical lattice, cos(2:) — + 1, the Hamiltonian ([7]) with 
the detuning pT|) is identical to the Landau-Zener model 
[2^ 1. A fast sweep (large Aq) implies non-adiabatic evo- 
lution and a large change in AE{t). Another possible 
reason for breakdown of adiabaticity is the non-linearity 
of the Gross-Pitaevskii equation. It is known that strong 
non-linearity (large values of g) tends to decrease adi- 
abaticity ^37i] . On the other hand, for certain systems 
moderate non-linearity can even increase adiabaticity 
[3^ . In Fig. [8] we present examples of the relative en- 
ergy, both for different sweep velocities Aq and different 
strengths g. For the solid/dotted/dashed lines we use 
zero scattering g = 0, but change the sweep velocity 



IV. CONCLUSION 

We have presented a thorough numerical analysis of 
spinor BEG driven through a critical point. The pres- 
ence of a first order FT for the interacting system was 
motivated by first studying the spectrum of the single 
particle Hamiltonian. Since the spectrum is non-gapped, 
the model could be sensitive to non-adiabatic excitations. 
Despite this fact, it was shown by solving the full Gross- 
Pitaevskii equation that the FT is still clearly visible even 
at relatively fast sweeps through the critical point. The 
effects of atom-atom collisions (Fig. 2]) as well as har- 
monic trapping (Fig. [7]) were also analyzed. Both were 
believed to smear out signatures of the FT. However, 
our findings demonstrate that the FT is resilient to both. 
Due to limitation of computational power, most results 
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